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Qv^ , In this paper we derive the exact expression of the four-spinon 

contribution to the dynamical correlation function of the spin 5" = 
1/2 anisotropic (XXZ) Heisenberg model in the antiferromagnetic 
^^. regime. We extensively study its isotropic {XXX) limit and derive 

r~| , perturbatively the Ising one. Our method relies on the quantum 

affine symmetry of the model, which allows for a systematic diago- 

k> ■ nalization of the Hamiltonian in the thermodynamic limit and for 

^ • an exact calcualtion of matrix elements of local spin operators. In 



fact, we argue that the familiar criticism of this method related to 
the complication of these matrix elements is not justified. First, 
we give, in the form of contour integrals, an exact expression for 
the n-spinon contribution. After we compile recently found results 
concerning the two-spinon contribution, we specialize the n-spinon 
formula to the new case n = 4. Then we give an explicit series rep- 
resentation of this contribution in the isotropic limit. Finally, after 
we show that this representation is free of divergences, we discuss 
the Ising limit in which a simple expression is found up to first order 
in the anisotropy parameter. 



Introduction 

Because of their well pronounced strong quantum betiavior and tiigtily nontrivial 
many-body interactions investigated in neutron scattering experiments on ferromagnetic 
compounds, and because of their rich mathematical structures, quantum spin chains have 
been the subject of intensive study over the past seven decades [1-19]. In this respect, one 
important quantity that plays a crucial role in these magnetic neutron scattering experi- 
ments is the dynamical correlation function (DCF) of local spin operators. This is because 
the cross section of such a scattering process is directly proportional to this DCF. In the 
past, most of the theoretical work related to such experiments has relied on approximation 
schemes [20-22]. 

However, certain one-dimensional compounds can be adequately described by models 
like the anisotropic Heisenberg (XXZ) chain, which has been shown to be amenable to 
exact manipulations [7,14]. Indeed, many static physical quantities of the XXZ chain have 
been determined exactly. But in contrast, dynamical physical quantities, and in particular 
the DCF, have remained elusive, due mainly to the complicated structure of the Hilbert 
space, which is built on " spinous" with nontrivial statistics and interactions. This is partic- 
ularly the case in the antiferromagnetic regime and the thermodynamic limit. Now as we 
understand it, the main reason why such a model is tractable exactly is its invariance under 
the quantum afRne symmetry. The recognition and systematic exploitation of this symme- 
try has led to the exact diagonalization in the thermodynamic limit of the corresponding 
Hamiltonian [23]. The Hilbert space is well defined and consists of multi-spinon particles. 
Moreover, using the technique of bosonization of the corresponding quantum affine algebra 
[24,25,26], one is able to calculate exactly matrix elements of local spin operators. This is 
explicitly carried out by Jimbo and Miwa in [27]. 

Nevertheless, one main criticism plaguing such calculations is that the matrix elements 
of local operators are typically expressed in the form of complicated contour integrals which 
are hardly useful for practical numerical applications related to experimental data and/or 
exact finite-chain calculations. In fact, this is precisely the reason why most of the very 
recent developments have been restricted to the derivation, using these matrix elements, of 



only the two-spinon contribution to the transverse DCF of the XX Z chain and its XXX 
and Ising hmits: the comphcated contour integrals are simply not present [28-30]. Note 
however that in [31], the two-spinon longitudinal DCF, where two such contour integrals 
arise, has been investigated perturbatively around the Ising limit. As shown in Ref. [29], 
in the XXX limit, the two-spinon contribution accounts for much of but not all the total 
DCF. Therefore, it is natural to try to look into the next contribution, namely that of four 
spinous. As far as we know, no numerical or analytic result, even approximate, related to 
this latter is yet available in the literature. The main reason is that in any approximation 
scheme, it is exremely difficult to single out from the rest the contribution of just the four 
spinous, and this even for a finite chain. In fact, to the best of our knowledge, the first 
attempt to tackle in a systematic way the case n > 2 is ref. [32] , where a formal expression 
for the exact n-spinon contribution to the DCF of the XXX chain is presented. But the 
contour integrals that arise there are not investigated thoroughly. 

It is the purpose of this paper to determine the exact expression of the four-spinon 
contribution to the DCF of the XX Z model in the antiferromagnetic regime and in the 
thermodynamic limit, and to discuss its isotropic XXX and Ising limits. In this contri- 
bution, a single contour integral arises with an infinite number of simple poles, integral 
that we express in the form of an infinite series. Even thus expressed, the four-spinon 
contribution in the XX Z model is still quite complicated, whereas its isotropic and Ising 
limits simplify considerably. Because of this, we concentrate our study on the XXX and 
Ising cases. In the XXX case, we show that our expression is free of potential infinities, 
and thus is well defined. In the Ising case, we derive a simple expression for the four- 
spinon contribution up to first order in the anisotropy parameter. We then compare it 
to the two-spinon contribution and discuss some discrepancies with results found in the 
literature. 

The paper is organized as follows. In the next section, following [23,27], we review the 
diagonalization in the thermodynamic limit of the XX Z model in the antiferromagnetic 
regime and, in particular, briefiy explain how it allows for an exact calcualtion of the matrix 
elements of local spin operators. In section 2, we define the DCF and use the completeness 
relation to write it as a sum over n-spinon contributions, n even, for which we derive a 



general coutour representation. In section 3, For the sake of completeness, we specialize to 
the case n = 2 and rederive all known results for the two-spinon contribution and update 
some of them. In section 4, we carry out a systematic study of the case n = 4 in the 
isotropic limit and discuss the Ising one. Finally, we devote section 5 to our conclusions. 

1. The anisotropic Heisenberg model in the antiferromagnetic regime 

In this section, we briefly review the diagonalization in the thermodynamic limit of 
the anisotropic Heisenberg XX Z model in the antiferromagnetic regime. We follow Refs. 
[23,27] where the quantum afline symmetry present in the model is exploited. This model 
is defined by the Hamiltonian: 

^ oo 

Hxxz = "2 E («+i + «+i + ^«+i) ' (1-1) 

n= — oo 

where A = {q + q~^)/2 is the anisotropy parameter. We consider the model in the anti- 
ferromagnetic regime characterized by A < —1 or equivalently by — 1 < g < 0. Here cr^'^'^ 
are the usual Pauli matrices acting non-trivially at the n^^ position of the formal infinite 
tensor product 

W = ---V®V®V--- , (1.2) 

where V is the two-dimensional representation of Uq{sl{2)) quantum group, with basis 
{v^,V-}. The main point of Refs. [23,27] is that the action of Hxxz on W is not well 
defined due to the appearance of divergences. However, because this model is symmetric 
under the quantum group Uq{sl{2)), its Hilbert space, which is free of divergences, is 
identified with the following level Uq{sl{2)) module: 

J^ = Y,ViA^)0ViA,r. (1.3) 

Here A^ and V{Ai);i = 0,1 are level 1 t/q(s/(2))-highest weights and t/5(s/(2))-highest 
weight modules, respectively, whereas V{Ai)* are dual modules defined from V{Ai) through 
the antipode. The module V{Ai) is identified with the subspace of the following formal 
semi-infinite tensor product of V^s: 

X = ---V®V®V, (1.4) 

4 



and which consists of all linear combinations of spin configurations with fixed boundary 
conditions such that the eigenvalues of a^ are (—1)*"'""^ in the limit n ^ — oo. Similarly, the 
module V{Ai)* is associated with the right semi-infinite tensor product of the V^s. What is 
impressive about this method is that it allows for a Fock space (or spinon-particule picture) 
representation of the Hilbert space, that is, this space can be created from a vacuum state 
through the successive actions of creation operators. This is achieved in terms of what is 
called type II vertex operators \l/^~*'*(^) and which interwine the Uq{sl{2)) modules as: 

vl/i-^'^(0 : V{A,) ^ V{0 ® V{Ai-,), (1.5) 

where V{^) is the affinization of the two-dimensional representation V in such a way that 
it is isomorphic to F §>> C[^, ^~^], with C[^, ^~^] being the algebra of polynomials in ^ and 
^~^. Here, ^ is a spectral parameter useful in parameterizing the energies and momenta 
of the spinous, which are by definition the eigenstates of the Hamiltonian. This spectral 
parameter lies on the unit circle. More specifically, the creation and annihilation operators 
are constructed as follows: decompose \l/^~*'*(^) as 

£ = ±1 

and let (\l/^~*''(^))* be the conjugate operator to \l/^~*'*(^). Then the following important 
commutation relations have been shown in Ref. [23]: 

(1.7) 

and which are the usual definitions of annihilation and creation operators respectively. 
Moreover, the vacuum state in the sector i is defined to be 

|0 >,= c-'/'(-g)^/2-^P„ (1.8) 

where Pi is the projector of y(Ao) ©F(Ai) on V{Ai), p is the grading element of Uq{sl{2)), 
and ci = {q'^ ; q"^) oo (see relation (1.13) for the definition of {y;x)oo)- The multi-spinon 
eigenstates are then systematically constructed as 

l^n, ■ ■ ■ 6 >.i ,....„;.= Cr^'K^ (^n) ■ ■ ■ K (^l) |0 >^ • (1-9) 

5 



Here C2 = (?^; ?^)oo/(?^; Q'^)oo and for clarity we have omitted the indices refering to the 
sectors on which the \1/* (^j)'s are acting, and the e's are the projections of the ;2-coniponent 
of the total spin. The completeness relation in JF reads then [27] 

I-2Z2Z ZZ ;^/ 2^'"2;^'^"''"'^'^'"'-''^'^ '''^'■••''"^^''■■■'^"'■ 

(1.10) 
The actions of Hxxz and the translation operator T (which shifts the spin chain by one 
site) on the spinon states are given by: 



T\Cir ■ ■ ,Cn >^ = l[T{C^)-^\Cl,■ ■ ■ ,Cn >l-^ , T|0 >,= |0 >i_, ; 
i=l 
n 



(1.11) 



where: 






<0 = r \'',' = e-'^(") , p{a) = ami— a) + ^/2 ; 



eiO = ^e:^logr(a = ^^ ^"'^^ Vl-fc^cos^(p) . 



(1.12) 



2q ^d^ ° '^' n 

In the above equations, e and p are the energy and the momentum of the spinon respec- 
tively, am{x) is the usual elliptic amplitude function with the complete elliptic integrals 
K and K', and moduli k and k'. Also: 

q = —exp(—TvK'/K) ; 
C = ^e'" ; 

00 

{y;x)^= Ylii-yx"") . 

n=0 

Furthermore, cr^'^'^(t) at time t and position n is related to a"o'^'^(0) at time i = and 
position n = via: 

<'^'^(t) = exp(itifxxz)T- V;^'^'^(0)T'^ exp(-itifxxz) • (1.14) 

At this stage, one should have expected the existence of "type I" vertex operators. Indeed, 
the latter denoted by $^~*'*(^) are also interwiners of the same modules as those of type II 
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vertex operators but in opposite order. Let us remind that the order of modules appearing 
in tensor produtcs is relevant in non-cocommutative algebras such as Uq{sl{2)). More 
precisely, they intertwine these Uq{sl{2)) modules as: 

$i-'''(0 : F(A.) ^ F(Ai-.) ® V{0. (1.15) 

They are similarly decomposed as 

$i-M(^)= J2 $i-^'^(0®t;,. (1.16) 

e=±l 

"La raison d'etre" of such type I vertex operators is that they allow for the translation of 
the action of local operators from that on the space W to one on JF. To be specific, any 
local operator can be written in terms of two-dimensional unit matrices E^-^^^, which in 
turn can be represented in JF as: 

^MiM. =C2$;,(l)$M2(l)®id. (1.17) 

Using the properties of the scalar products defined on t/q(s/(2))-highest weight modules 
it has been shown in Ref. [23] that the matrix elements of local operators can be then 
expressed as traces over highest weight modules of products of type I and type II vertex 
operators. For example the matrix element of the unit local operator E^^^^ acting on site 
1 is equal to 

^.|p I, , ^ _ l-n/2 ^y (A. (9-^^+-^;, (1)^* , (1)^L {in) ' ' ' ^t, (6)) 

I ^ U|£/^^^2lSn, ■•■,?1 >£„,-••, ei;i ^2 / 2p+A 

(1.18) 

Unfortunately, these traces are over infinite-dimensional modules and hence cannot be 
easily evaluated. What is so special about the anisotropic Heisenberg model in the antifer- 
romagnetic regime is that the quantum affine algebra Uq{sl{2)), its highest weight modules 
y(Aj), and both its type I and type II vertex operators can be bosonized in terms of a 
set of simple harmonic oscillators satisfying the usual Heisenberg algebra*. Therefore, this 



* Heisenberg not only defined the model, he also defined the right tool for its exact 
resolution, though he did not realize it at that time 



problem is translated to a familiar quantum field theory problem of computing the Green 
functions of a free bosonic field. Using the Wick theorem and other procedures the above 
traces have been explicitly computed by Jimbo and Miwa [27]. From this brief summary 
it is tempting to postulate that whether a model is explicitly integrable or not depends on 
whether it is bosonizable or not. 

2. The n-spinon transverse DCF of the XX Z model 

First, we define one of the two (equal) transverse components of the DCF in the case 
of the XXZ model. The reason why we focus only on the transverse DCF is that it 
involves the least number of contour integrals and, in particular, the transverse two-spinon 
one does not involve any. This is to be contrasted with the longitudinal DCF which, even 
in the simplest case of two spinous, involves two complicated contour integrals and this 
number increases with the number of spinous. We define this component as: 

/oo 
dt J2 e^^^'+'^'^h < 0|a+(t)(7o-(0)|0 >, , (2.1) 

-°° mez 

where w and k are the neutron energy and momentum transfer respectively, and i corre- 
sponds to the boundary condition. Using the completeness relation, the n-spinon contri- 
bution is given by: 

^'- ^z.hJ ^^'^^ ^^<- k (2-2) 

X i^-a < 0|cr(|(0)|^n, •••,6 >e„,...,ei;i+n i;ei,...,e„ < ^1, ■ ■ ■ , ^nWo {0)\0 >i . 

It can be rewritten in the following more tractable way: 

ei, ••-,£« meZ j = l 

X (^, <0|a+(0)|^„,...,a >e„,...,ei;z ^;el,...,e„ < 6, • • • , ^nl^o" (0) |0 >, 

+ ei^+T^U ^«^)) i_, < 0|ao+(0)|^^, . . . , 6 >e„,...,e,;l-. 

z;ei,...,e„ < 6, • • • , ^n^o" (0) |0 >, ) . 

(2.3) 



The non- vanishing matrix elements have been computed in [27]. In particular, they satisfy 
the following relations: 

.<0|ao"(0)|^n, . . . , a >.„,.. .,ei; ^ = l-^< 0|(T+(0)|a, . . . , ^ >-.„,.. .,-ei; l-^ 

= ,< 0|ao+(0)| - q^^\ . . . , -qC' >-.!,...,-.„; ^ ; (2-4) 

i;ei,...,e„<6,---,'^n|f7(7(0)|0>i = i< 0|(J(7(0)| - q^l, . . . , -q^r, >-£i ,...,e„; i • 

We ease slightly the notation by writing: 



^l,...,..(en,...,6)=^<0|ao+(0)|a,..., a >.„,...,.,; 



(2.5) 



and we have: 



X 



l,...,e,(en,...,a)=(5^"^^,^,_J- 



ni<j<i'<n7K'/%) 



nLi(-?'"i;?^)oo(-?^'U- ;g 



3., -1.^4^ 



n 



j-(l+e,)/2-i 
3 



n£^n(''«7'-«-v))n(«.-'-"7')) 

1 



{uj/vf,q^)oo{q '^vi/uj]q'^ 



W{-q ^vi]q^)^{-q^/vi]q' 



3,i 

\Y V£'{vi' -q~'^V(){v(,/v(,']q^)^[q^vii/v(,;q^)ooG^{v,u) 
In the above equation, we have used the following definitions: 



(2.6) 



^3 = -^3 



Gti 



v,u 



^2 (3 — i)n 



2-n+4:i [\jyi_ 



C+,=5,,^-l{-q)-'^^'^-\l-q-')--\q';q'U{q';qr^-'p-{q); (2.7) 



7(0 



1. ^4. „4\ 



(q ^9 ;9 )oo(^ ;r;9 



(g^;g^;Q^)oo 
p{q) - 



(g^;g^;g4)oo 

Moreover, t is the number of elements in the set L, which is defined by: 

n 

L = {j, s.t. ej = +1, and /J^i = —2, 1 < j < n} . 



(2.^ 



Finally, in (2.7), each contour Ci includes the poles q Uj, n > only, excluding the poles 
q^-^n^,^ n > 0. 

Now, restricting ourselves to the first Brillouin zone (i.e., < /c < 27r), integrating the 
delta functions, and keeping track of the Jacobian factors, we find: 

S^^ [w, k-n) = Cn f 7--— . . . 7--— > ^ , 2.9 

/ 27rz^3 2m^ri,'^, J{Pi,P2) 

where we have: 

x^:„,...,e,(c\---,er'); 

(2.10) 



J{Pi^P2) =1 sin(2pi)A/l - /c2cos2(p2) - sin(2p2) a/I - k'^ cos2(pi)| ; 

a. 



4nl{q - q-^)k'^K^ 

In equation (2.9), the sum over the pairs (^1,^2) is a sum over all the solutions to the 
energy-momentum conservation laws: 

n 

«; = e(a)+e(6)+Xl^(^^) ; 

^"^ (2.11) 

k^-p{Ci)-p{C2)-Y.P(^j)- 

J=3 

The physical DCF is obtained by summing over both sectors i = 0, 1. Finally, let us 
mention that the isotropic limit of this integral representation of the n-spinon contribution 
has been presented in Ref. [32]. 

3. The two-spinon DCF and its isotropic and Ising limits 

In this section, we discuss the case n = 2. This is particularly interesting because 
no contour integral arises, that is, the set L in equation (2.8) above is empty. We briefiy 
rederive the known results [28-30] and, more importantly, comment on some new properties 
of theirs. The two-spinon DCF takes the following simple form: 

Sl^+-{w,k-n)=C2^^jj^^. (3.1) 

10 



Note that 5*2 vanishes outside the domain of 2-spinon continuum set of states. The rest of 
the notations is as follows: 

o,^{qei)OAQQ 



7±(0 



I oo 



(^g4(_^)3±l. ^4. g4)^(^-l(_g)3±l. ^4. ^4) 



J{Pi-,P2) = I sin (2^1)^1 - /c^cos^ (^2) - sin (2^2) yl - /c^ cos^ (pi)| ; 
^2 



C2 = 



(3.2) 



2{q-q-^)k^K^ ' 
The measurable DCF is obtained by summing over both sectors, i.e.,: 

S+'iw, k-n) = S^'+-{w, k-n) + S^'+-{w, k-n) . (3.3) 

Finally, for fixed w and k, (ai,a2) is a solution to the energy-momentum conservation 

laws: 

w = e(ai) + e(a2) ; 

(3.4) 
k = -pi - p2 ■ 

3. a The isotropic limit 

As it is clear from the two equations (3.1) and (3.2), the 2-spinon contribution in 
the anisotropic model is still quite untractable for practical purposes. However, certain 
physical limits of that result are quite simple, in comparison. Here we briefly rederive the 
known result of the isotropic limit [28]. Let p = —q be the nome of the various theta 
functions that appear in 5*2 ~. The isotropic (XXX) limit is obtained by first making 
the redefinitions: 

(3.6) 

and letting e — > 0+. In this isotropic limit, the real number (3 is the appropriate spectral 
parameter. Also, in this limit, we have the following asymptotic behaviors: (see also Ref. 

11 



[27] 






7±(-r) 



2^ (/;/)oo(l-p^)^/^^^/^A±(/3) . (3.7) 



r(l/2±l/4)A±(z7r/2) 



(-P«^;p^)„(-prV^^ - (i-^'^x^^^''^)^ 



/oo ~ 



r(l/2-/3/7rz)r(l/2 + /3/7r«) 



where: 



. ,. I p^^ sinh^t(l--^)exp(Tt) \ 

A+f^j = exp - / dt — — ^ — -— 

^^ ' ^ ' Jq tsmh(2t)cosh(t) / 



(3i 



I 7 

|A±(^j| -exp|^ y^ dt tsinh(2t)cosh(t) J ' 

with z = X + iy and x, y being real. Using the fact the for generic k, the two pairs 
(P17P2) and {p2,Pi) satisfying the energy-momentum conservation relations have equal 
contributions to S2 , then this latter reduces to the following expression [28]: 

,.-,„,_„_, 1^,(^^(P^^^_ (3,8) 

with: 

w = -7r(sin(pi) + sin(p2)) ; 

k = -pi - p2 ; 

(3.10) 
Wu = 27rsin(A;/2) ; 

Wi = T{\ sin/c| . 
Let us mention that 5"^^ (w, k — ix) has been expressed in [29] in terms of w and k through 
the following relation: 



9 2 

wt — wf 



/9bl)-/3b2)=2cosh-^/^ L . (3.11) 

y ty^ — wf 

Also, recall that S2 ~ vanishes outside the domain limited by Wu and wg. Summing up 
over the two sectors and over the two possible contributing pairs of spinous, and using the 
relation: 

a^ = ' , (3.12) 



There is misprint there which is corrected here. 
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we find: 

S2^{w,k — n) = 8S2 ~{w,k — Tx), ii = x,y,z. (3.13) 

Note here that we have included the longitudinal component of the two-spinon DCF with 
the two transverse ones because they are all equal in this isotropic limit. For further 
applications of this result, see [29]. In section 4, we discuss in detail the four-spinon DCF 
in the isotropic limit. 

3.b The Ising limit 

The other practically interesting limit of the anisotropic DCF is the Ising limit, which 
is already known in the literature [33]. Here we briefly rederive it in the context of the 
present theory, see also [30], and compare it to the perturbative result of Ishimura and 
Shiba [33]. 

The Ising limit is the limit in which the interaction term ^ ^^ cr^a^j^i in the Hamil- 
tonian becomes dominant. It is therefore obtained by first rescaling the Hamiltonian as 
-falsing = Hxxz /\^\^ and then letting the no me — g = e ^ 0+. This means the energy 
e of a spinon and hence the factor J{pi,P2) in eq (3.1) above must also be rescaled as 
eising = e/|A| and Jising{pijP2) = J{PijP2)/\^\, respectively. Using also the following first 
order expansions in the limit e ^ 0: 

K~^(l+4e + 0(e2)); 

2 (3.14) 

A;~4ei/2(l-4e + 0(e2)) ; 
we find, for fx = x,y 

Sifngiw, k-Tx) = 250'+-(«;, A; - 7r) + 23^^+' {w, k - it) ■ 

_ V4|yp-(«;-2)V {w-2){l + cosHk)) . (3.15) 

|y|2 I' 4ecosi/cj 2cos2(A;) J ' 



with: 



e{p) ~ 1 — 4ecos(2p) ; 

w = e{pi) + e{p2) ~ 2 — 8ecos(A;) cos(/c + 2pi] 
V = 4ecos(/c) ; 
\w-2\ <2\V\ . 
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(3.16) 



Note that, except for the third one, aU other factors in the above formula agree with those 
of Ishimura and Shiba [33]. These authors have used a completely different method, i.e., 
the Green function technique. Also, note that both formulas completely coincide in the 
particular cases of w = 2 and k = n. 

The most important difference is that in our expression, the order of the pole at 
k = 7r/2 is increased by two at the zeroth order in q. It is our opinion that the disagreement 
between the two formulas is due to the fact that the authors of [33] have used a dispersion 
relation only up to first order in e. However, following Ref. [21], the DCF can be written 
as: 

S+~{w, k-7v)= D{w, k-Ti) M{w, k - 7r)/2 , (3.17) 

where we have: 

1 

D{w, k-n) = 2\A\/\dw/dpi\ 



8|sin(2pi + k)cos{k) e-2sin(4pi + 2k) cos{2k) e^l 



1 , , cos(2iDi + k) cos(2/c) , 
(1 + 4 ^—^ 7T-. ^ — - e) 



8e| sin(2pi + k) cos{k) \ cos(/c) 

(3.18) 

From the above relation, it is clear that a second order dispersion relation in e contributes to 

zeroth order in the DCF. Therefore, it seems that this second order term in the dispersion 

relation is missing in the Ishimura- Shiba formula whereas it is present in our expression. 

4. The four-spinon contribution to the DCF 

Let us now turn to the main subject of this paper, namely, the four-spinon contribu- 
tion. As mentioned in the introduction, as far as we know, no previous study of this case, 
whether numerical or analytic, is available in the literature. After we process a little more 
the anisotropic case, we quickly focus our study on the isotropic case, and then we discuss 
the Ising limit. As the reader will rapidly realize, the reason is that in the anisotropic case, 
the expression of the four-spinon contribution, while simplifying somewhat, is still difficult 
to manipulate. 

We let n = 4 in (2.3) and in the subsequent relations of section 2. The set L has one 
element and so, there is one contour integral to perform. The four-spinon contribution (that 
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we denote 5*4) is defined in terms of the matrix elements X|(^4, . . . , ^1), where £ = 1, . . . , 4 
corresponds to the position of the + sign in the sequence {ei, £2, £3, £4} that enters in the 
definition of the set L, equation (2.8). This matrix element can be written as: 

^i ni<j<fc<4 7(^fc/'Uj) -p-,- j-(l+^^)/2-^ 



]lj=i{-'iUj] q^)oo{-q^Uj; q 



4:\ 

00 j^l 



where: 



dV£ T-r Uj — V£ T-r qUj — q ^V£ 



U^-^U 



2niv£ -'■-'- UiV£ -'■-'■ UnV£ 



4 
X 
J 



4 1 _ 



(4.2) 



The contour is as described in section 3, after equation (2.8) and G~^''^{v£) simplifies to: 

G+'\v£) = {-v£ye,s{-q-'+^^ ^' )c+, 

U1U2USU4 ' (4.3) 

ct = (-?)'-''(! - q-')iq'; q')oo{q'; /)Lp'(?), 

Because all the poles are simple, it is therefore possible to perform the contour integral 
present in the above expression of l£. Being careful about the zeroes that annihilate some 
of the poles, we find the following series expansion for it: 

4 



j=l m=e{j-e) i<i i>i 



'-U-'] 



"" n (^., 4)^(^-2+4^^. ^4)^ nr^^(^4. _ ^) (4-4) 



X 



(_g-l+4m„.;Q2)^(_^^-1.^2)^^+,^(^4m^^.)J-[2^-1^^4^_l_2.^^-l. 



(?^?4)oo(?-2+4-;Q4)^J|-^(g4._l) 

The recipe to derive a series expansion for 5*4 is to substitute this expression for l£ 
back into the expression of Xj, and then substitute the latter in the expression of 5*4. The 
general term we thus find is quite complicated and therefore we prefer to omit reporting 
it in this article. We rather focus our attention on the isotropic limit (and later on we 
discuss the Ising limit) where major simplifications occur. 

4.1 The isotropic limit 
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The isotropic limit is performed by first redefining Ui = —^f = e^'^^*'*'^ and p = —q = 
e~^ and taking tfie limit e ^ 0+. Using the following relations, which give the leading 
terms in this limit of the different expressions involved in this calculation: 

p'^iQ) = 7+(9~^) -^ liu) = 7_('u) ; 
u = —4 = e '^^ ; 

ui - U2q^'^ ~ — (/3i - /32 + 2mm) ; 
in 

ui - U2q-^+^'^ ~ -(A - /?2 + 7rz(2m - 1)) ; 

^7r (4.5) 

/ 4x / 3 -1 4x (9^;9'^)L(l-9^) 

{-qu]q )oo{-q u ;q )oo -■ 



■^ V4 271-i/-^ V4 ^ 271-i/ 
(_g-l+4-^; Q2)^(_g3-4^^-l. ^2)^ ^ _(^2. ^2)2^(1 _ ^2) ^^^j^^^^ . 

we find after a lengthy algerba that the isotropic (XXX) limit of 84^ simplifies to: 
a rather manageable expression. The notation is as follows: 



JP 



(4.7) 



/(/3i,/?2,/?3,/34)= n i^-(/5/-/3. 

l<i<j'<4 

r ^ I . 

^ 3x29r(l/4)8|A_(z7r/2)|8 ' 
VF = zw + 7r(sinp3 + sinp4) ; 

V7„ = 27r|sin(K/2)| ; 

K = k+p3+p4 ; 

cot pj = sinh (3j , — tt < Pj < , 

and gr^ is to be defined shortly. Also, for fixed W and K, the sum ^(p^ p^) ^^ °^^^ ^^^ ^^° 
pairs (pi,P2) and (p2,Pi) solutions to the energy-momentum conservation laws: 

W = — 7r(sinpi + sinp2) ; 

(4.8) 
K = -pi - P2 , 
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and which are in fact given by 

N f K , W . K / W .\ ,, , 

{PuP2) = [- — + arccos (- — — ^) , -— - arccos (- — — ^) j . (4.9) 



The rest of this subsection is devoted to the discussion of some of the properties of 
the integrand of 5*4. We start by discussing the function (7^, which is given by: 



g, = (-1) (2.) 5:cosh(2.p,) J2 U..,n-'sinHnp,^) 

^ -A- r(m - i + ipj^) 



(4.10) 



-^ r(m + 1 + ipji) 
where we have set (3j/2tt = pj and pj — pi = pji. It can be rewritten as: 



4 



1 -p-|- A/Sinpismpj ^-^ -p-j- 1 



j=i ^^^Po ^^j sin \^^) ^^eo-z) */^ 

r(m - I + ipji) 



X 



n 



A r(m + 1 + ipji, 

expression in which we have used the foUowing: 

cosh(27rpj) sinpj = —1 ; 
sin (^^) 



(4.11) 



sinh(7rp,,) ... 

^smpiSinpj (4.12) 

cosh(7rpji) = — 



sin(^) 



y/sinpi sinpj 

Both expressions for (7^ are useful in the sequel. 

We first discuss the behavior of gn in the region where two /Sj's (equivalently two Pi's) 
are equal. We find that gi is finite for all I. The same is true in the regions where three 
of the /3i's or all the four are equal. Note that for a given (7^, each term in the X]7=i ^^^ 
be divergent, but when the divergent pieces are put together, they cancel one another. To 
be specific, consider for example g\ in the region p\2 very small (the other pji being kept 
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finite). Then we have: 



gi/[2Tv] = cosh(27rpi) > rmr . , , , . , , r f : r 

m=l 1=1 

1 u/o ^ V^ 2("^ + ^P23)(m + zp24) Ar(m- i + ^p2^) (4-13) 

H cosh(27rp2) / ^t^tt — — — -— 

P21 ^^ sinh(7rp23) sinh(7rp24) fJi r(m + 1 + ip2i) 

+ regular . 

Now replacing p2 by pi in any regular part and using pi2 = — P21, we find that the singular 
pieces cancel each other and the resulting nonvanishing function is regular. The same is 
true in all the other cases.* 

We now turn our attention to the asymptotic behavior of agrn, the general term of the 
series in gg. The behavior of agrn at large m can be obtained if we use the asymptotics of 
the F-function. We find that at large m: 

1 1 , — 11 

T{m- - + ipji) ~ r(m--) ~ V27rexp [(m - 1) In (m - -) - m + -] ; 

T{m + l + ipji) ~ r(m + 1) ~ v^exp [(m + -) In (m + 1) - m - 1] ; (4.14) 



-I— r 1 



2' 
3 



We gather all the pieces together and we arrive at the following asymptotic behavior: 






m^ 



which leads to lim^^oo |«m+i/«m| = 1~- The series is thus convergent for generic pj, 
though not 'very' rapidly. 

Finally, we discuss the behavior of the integrand of 5*4 when one of the four p^'s goes 
to ±00 ( equivalent ly, pi -^ — tt or 0). To be specific, and because in ^'4 there is an explicit 
integration over ^4, we consider explicitly the case p4 -^ +00 while keeping the others 
finite. Note that the conclusion we arrive at is the same for all the other momenta. We 



Since |A_| goes to zero in these regions, see [29], then the whole integrand of 84^ has a nice regular 
behavior there. 
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first discuss the behavior of gi and then we incorporate / into the picture and discuss the 
behavior of / Yli^i Ifi'^P- 

To start, let us consider for example (7i/(27r)'^. Explicitly, it is equal to: 



n,=2("^+^piO 






uL2i^ + ip2i 



+ cosh 27rp2 > rt 1 ■ x.( \ 11 r 

i-'i ^^/2 ^ sinh(7rp2z) fj^ r 



nJ=2("^+^P3* 



+ cosh 27rp3 > :pf . M - 



+ cosh(27rp4) ^ — 

m=l ^^ 



^^=2K + ^P4^ 



n 



m - 2 +ipii] 



m + 1 + zpii 



TO 



2 + ip2z 



m + 1 + zp2i) 



(4.16) 



m - 2 + «P3iJ 



m + 1 + zp3i) 
m - i + zp4i) 



TO + 1 + Zp4i 



;^_^4 TT 1 sinh(7rp4i) ^1 r( 

For illustration, denote the term containing cosh(27rpi) by the superscript (1). Using the 
asymptotic behavior of T{z) at large 1^1, we have: 

r(m - ^ + ipiA) 1 



r(m + 1 + ipi4) 
We then re- write this first term (1) as: 

cosh(27rpi) >^ (m + ipi2) (m + ipia) 



(m + 1 + ipi4 



(4.17) 



-p-j- r(m - i + ^pli) (m + zpi4) 



TT 1 sinh(7rpi4) ^^ tt ^ sinh(7rpi2)7r ^ sinh(7rpi3) jL"^ r(m + 1 + ipu) (^ + 1 + ipi^)i 

(4.18) 

The series above still converges as 1/m"^, thus keeping only a finite number of terms in it, 

we can neglect m in front of pi4 and we get the asymptotic behavior of this first term: 



9i 



(1) 



P4 ^ e-^^4 i? 



(4.19) 



(We have neglected pi in front of p4 and we have extracted the relevent exponential from 
the sinh term). Here the function i? is a generic notation for whatever is known to be 
regular. We work out the other three terms in the expression of gi/{2'K)'^ in a similar 
manner and we find that they behave asymptotically in the same way as this first term 
behaves (except the fourth one which goes to zero faster). Thus we can write: 



^1 ~ P4 ^ e-^P^ R . 



(4.20) 
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In fact, the same pattern happens for the other (7^'s and we arrive at: 

J2\9i\^ e'^""'* R- (4.21) 

£=1 P'^ 

Equivalently, the asymptotic behavior of gi can be derived in the p-variables, and in the 
hmit p4 -^ — 7r, we find: 

^ ln(-smp4) 

Now we put into the picture the contribution from /. Using the notation of [29], / is 
given by: 

/ = exp[-/(ti2) - /(tl3) - /(tl4) - /(t23) - /(t24) " I{h^)] , (4.23) 



where: 



-/(t) = C2/2 + In (t^/2 sinh — ) + h{t)/2 , (4.24) 



and tji = 4pji*. One checks exphcitly that: 

h{t) ^0 as t^ +00 , (4.25) 

which means that when p4 -^ +00, 

/ ~ pTe'^P* R , (4.26) 

where we neglect the other p's in front of ^4. Using the behavior (4.23) above, we find: 

4 
/ Y.\9i\^ ~ pTe'^'^'R. (4.27) 

Had we used the pi variables, we would have obtained the behavior: 



/ 1 1*1' ~ ./Mz£Ea) R . (4.28) 

^ y -smp4 



* Note that in Ref. [29], there are few misprints regarding the definition of I{t) that 
are corrected here. 
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The convergence of the integrand is best seen if we use (ip4 = —iTidp^j cosh (27rp4) to write: 

4 

dv^ f Yl \9e\^ ~ P^^^ ^""''' ^ "^P^ ■ (4.29) 

£=1 

It is interesting to note that this asymptotic behavior of the integrand of 5*4 resembles that 
of 5*2 in the vicinity of the lower boundary w = wg, and which also corresponds to the 
limit where one of the two momenta goes to or —n. 

4.b The Ising limit 

We finish our study by considering the behavior of 5*4 in the Ising limit. This limit is 
performed as in the 2-spinon case, namely by first rescaling the energies of the four spinous 
by |A|~^ and thus the Jacobian in the expression of S'4, and then letting q -^ 0~ . We 
return to equation (4.4) above where a series expansion of /] is given. Because of the factor 
Y[i<eiQ^J^ ^ (l'^^~^u~^), it is clear that the term with i = 4 and m = is the leading one 
in powers of q. Explicitly, we find: 

'-' "*J ' (4,30) 






1\ 



Furthermore, from the definition of a in terms of ^, see equation (1.5), we have the following 
relation between the physical momenta pi and the spectral parameters a^ (up to first order 
in q): 

TT 

tti =Pi -2Qsin(2pi) - - . (4.31) 

Now, using the following expansions: 

7-(^)7-(i) ~ 4sin2(ai-a2)+0(Q2 
1 + 0(q), 



-72''7-(,72 . 

, 2 2,4 (4-32) 
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We find that h^ + h^, which appears in the numerator of the integrand of 5*4, greatly 
simplifies to: 

/^O , ^1^2^ 2 '^ ( i\m+n-2 COSJPm - Pu) Smjprn - Pn) Sm{p^ + /c/2) sin(p„ + fc/2) 

X Yl sin(pj -pj,) , 

l<j<j'<4 

(4.33) 

where: 

rh = min{m, n) ; 

fi = max{m, n) ; 

(4.34) 
m = min{{l, 2, 3, 4}/{m, n}) ; 

n = max({l, 2, 3, 4}/{rM, n}) . 

For example if (m, n) = (1, 3) then 1 = 1, 3 = 3, m = 2, and n = 4. Moreover, using the 
following expansion of the denominator of the integrand of S^: 

C 

— r^ 

siiv{2pi)Jl - k'^cos^{j)2) - sin(2p2)\/l - k'^cos'^{pi) , , ^,.^ 

? (1 _ ^^ cos(pi-p2) ^ ^ 

2q cos(pi + P2) sin(pi - P2) cos(pi + P2) 

we arrive finally at the following simple integral for the leading term in g of 5*4: 

96^-^;^ ixm+n-2 /"° /"" , , sin(p^ - p^) sin(p^ + /c/2) sin(p„ + /c/2) 

-34 ~ ^ > (-1J / / dpsdp4^^^- -— — r ■ -— 

'^^^i J-nJ-n sm{p^-pf,)\sm{pi-p2)cos{pi+p2)\ 

XCOs{pm-Pn) JJ Sm{pj - Pj') . 

l<j<j'<4 

(4.36) 
Let us recall that pi and p2 are given in terms of pa and p^ through the energy-momentum 
conservation laws. Consequently, the ratio of the contribution of four spinous to that of 
two spinous is of the order of g^, thus confirming all previous numerical studies which show 
the dominance of the contribution of two spinous. 



5. Conclusion 
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In this paper, we have studied the exact spinon contributions to the dynamical corre- 
lation function (DCF) of the anisotropic XX Z Heisenberg model in the antiferromagnetic 
regime, and their isotropic and Ising limits. 

First we have given a general contour-integral representation for the n-spinon contri- 
bution, which we have used to rederive most of the presently known exact results on the 
two-spinon contribution 5*2 in the XX Z case and its XXX and Ising limits. 

Then we have focused our attention on the first nontrivial case for which no results 
(even numerical ones related to finite-chain calculations) are available yet, namely, the 
four-spinon contribution S4. In this context, we have shown that the contour integrals 
can be performed, yielding mathematically well defined expressions, free from divergences. 
This important result demonstrates that the usual criticism against the methods based on 
the quantum affine symmetry is unjustified. In fact, our results can be thought of as yet 
another instance of the usefulness of the quantum affine symmetry of the XX Z model. 

Finally, we touched upon the Ising limit of S4. We have shown that its leading contri- 
bution is of the order of q and thus, in this limit, the ratio of the four-spinon contribution 
to that of the two-spinon one is of the order of q^. It seems that the expression of 5*4, up to 
this order in the Ising limit, is the simplest result possible for the DCF beyond two-spinon, 
and thus we think it deserves further attention because it might provide new insights on 
the behavior of the highly nontrivial quantum mechanical nature of the Heisenberg model 
that is partly encoded in 5*4. Moreover, in this isotropic limit, we have found that though 
the series defining 5*4 converges, it does so quite slowly. Therefore, we think further nu- 
merical analysis is required to see what is the best criterion of truncation so that we get a 
satisfactory numerical evaluation, in view of any potential application to new experiments 
where the main focus is to be on the explicit detection and measure of any signature from 
four spinous. Indeed, as far as we know, no experiment thus far has effectively focused on 
this contribution. The main reason is that up to now, much of the attention of the theoret- 
ical work has been limited to the two-spinon sector, and even in this latter, our theoretical 
understanding has not been complete until we have started to exploit systematically the 
quantum affine symmetry of XX Z [23], thus allowing the derivation of the simple exact 
result of [28]. 
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From a general perspective, it would certainly be interesting to extend these methods 
to other models known to be invariant under quantum affine symmetry, or any other 
infinite-dimensional symmetry. This is already extensively implemented in conformal field 
theories. 
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